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Abstract. Zonotopal algebra is the study of a family of pairs of dual vector 
spaces of multivariate polynomials that can be associated with a list of vectors 
X. It connects objects from combinatorics, geometry, and approximation the- 
ory. The origin of zonotopal algebra is the pair (D(X), V(X)), where T>(X) 
denotes the Dahmen-Micchelli space that is spanned by the local pieces of the 
box spline and V(X) is a space spanned by products of linear forms. 

The first main result of this paper is the construction of a canonical basis 
for T>(X). We show that it is dual to the canonical basis for V(X) that is 
already known. 

The second main result of this paper is the construction of a new family of 
zonotopal spaces that is far more general than the ones that were recently stud- 
ied by Ardila-Postnikov, Holtz-Ron, Holtz-Ron-Xu, Li-Ron, and others. We 
call the underlying combinatorial structure of those spaces forward exchange 
matroid. A forward exchange matroid is an ordered matroid together with a 
subset of its set of bases that satisfies a weak version of the basis exchange 
axiom. 



1. Introduction 

A finite list of vectors X gives rise to a large number of objects in combinat- 
orics, algebraic and discrete geometry, commutative algebra, and approximation 
theory. Examples include matroids, hyperplane arrangements and zonotopes, fat 
point ideals, and box splines. In the 1980s, various authors in the approximation 
theory community started studying algebraic structures that capture information 
about splines {e.g. [Tl [12l ETj). One important example is the Dahmen-Micchelli 
space V(X), that is spanned by the local pieces of the box spline and their par- 
tial derivatives. See [551 Section 1.2] for a historic survey and the book [T3] for a 
treatment of polynomial spaces appearing in the theory of box splines. Related res- 
ults were obtained independently by authors interested in hyperplane arrangements 
(e.g. 

The space P{X) that is dual to T>(X) was introduced in [TJ[2T|. It is spanned by 
products of linear forms and it can be written as the Macaulay inverse system (or 
kernel) of an ideal generated by powers of linear forms [TT] . Ideals of this type and 
their inverse systems are also studied in the literature on fat point ideals [241125] , 
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In addition to the aforementioned pair of spaces (V(X), V(X)), Olga Holtz and 
Amos Ron introduced two more pairs of spaces with interesting combinatorial prop- 
erties |28| . They named the theory of those spaces Zonotopal Algebra. This name 
reflects the fact that there are various connections between zonotopal spaces and 
the lattice points in the zonotope defined by X if the list X is totally unimodular. 

Subsequently, those results were further generalised by Olga Holtz, Amos Ron, 
and Zhiqiang Xu [29] as well as Nan Li and Amos Ron [33]. Federico Ardila and 
Alex Postnikov studied generalised P-spaces and connections with power ideals [5] . 
Bernd Sturmfels and Zhiqiang Xu established a connection with Cox rings |41| . 
Further work on spaces of P-type includes (3) [5] [32j [42] . 

Zonotopal algebra is closely related to matroid theory: the Hilbert series of 
zonotopal spaces only depend on the matroid structure of the list X. 

It is known that there is a canonical way to construct bases for the spaces of 
P-type [H [2TJ HH1 [33 • The first of the two main results in this paper is an algorithm 
that constructs a basis for spaces of 2?-type. Two different algorithms are already 
known [9[ 115j. However, our algorithm has several advantages over the other two: 
it is canonical and it yields a basis that is dual to the known basis for the P-space. 
Here, canonical means that the basis that we obtain only depends on the order of 
the elements in the list X and not on any further choices. 

Our second main result is that far more general pairs of zonotopal spaces with 
nice properties can be constructed than the ones that were previously known. We 
define a new combinatorial structure called forward exchange matroid. A forward 
exchange matroid is an ordered matroid together with a subset of its set of bases 
that satisfies a weak version of the basis exchange axiom. This is the underlying 
structure of the generalised zonotopal 2?-spaces and P-spaces that we introduce. 

All objects mentioned so far are part of what we call the continuous theory. 
If the list X lies in a lattice (e.g. Z d ), an even wider spectrum of mathematical 
objects appears. We call this the discrete theory. Every object in the continuous 
theory has a discrete analogue: vector partition functions correspond to box splines 
and toric arrangements correspond to hyperplane arrangements. The local pieces of 
the vector partition function are quasi-polynomials that span the discrete Dahmen- 
Micchelli space DM(A). Both theories are nicely explained in the recent book by 
Corrado De Concini and Claudio Procesi [T7]. The combinatorics of the discrete 
case is captured by arithmetic matroids which were very recently introduced by 
Luca Moci and Michele D'Adderio 0135]. 

Vector partition functions arise for example in representation theory as Kostant 
partition function, when the list X is chosen to be the set of positive roots of a 
simple Lie algebra (e.g. [7]). There are also applications to the equivariant index 
theory of elliptic operators |18[ \W\ [20] . 

This paper deals only with the continuous theory. However, the two theories 
overlap if the list of vectors X is totally unimodular. We hope that the results in 
this paper can be transferred to the discrete case in the future. 

1.1. Notation. Our basic object of study is a list of vectors X = (x\, . . . , xjv) that 
span an r-dimensional real vector space U = M. r . The dual space U* is denoted by 
V. We slightly abuse notation by using the symbol C for sublists. For Y C X, 
X \ Y denotes the deletion of a sublist, i. e. (x±, X2) \ (xi) = (#2) even if X\ = Xi- 
The list X comes with a natural ordering: we say that xi < Xj if and only if i < j. 
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Note that X can be identified with a linear map I 1 *' — > U and after the choice 
of a basis with an (r x 7V)-matrix with real entries. 

We will consider families of pairs of dual spaces (V(X, -),V(X, ■)). The space 
T>(X, •) is contained in Sym(V), the symmetric algebra over V and V(X, •) is con- 
tained in Sym([7). The symmetric algebra is a base- free version of the ring of 
polynomials over a vector space. We fix a basis (si, . . . , s r ) for U and (ti, . . . , t r ) 
denotes the dual basis for V, i e. U(sj) = 5ij, where 6ij denotes the Kronecker 
delta. The choice of the basis determines isomorphisms Sym(£/) = R[si, . . . , s r ] 
and Sym(V) = R[ti, . . . ,t r ). For x G U, x° C y denotes the annihilator of i. e. 
i° := {/ E V : /(.t) = 0}. For more background on algebra, see Q~7] or |22j . 

As usual, xs : 5 1 — > {0, 1} denotes the indicator function of a set S. 

1.2. Matroids. An ordered matroid on iV elements is a pair (A,B) where A is 
an (ordered) list with N elements and B is a non-empty set of sublists of A that 
satisfies the following axiom: 

Let B, B' G B and be B\B'. Then, there exists b' G B' \ B s. t. (B \ b) U b' G B 

(1.1) 

B is called the se£ o/ frases of the matroid (A, B). One can easily show that all 
elements of B have the same cardinality. This number r is called the rank of the 
matroid (A, B). A set J C ^4 is called independent if it is a subset of a basis. The 
rank of Y C A is defined as the cardinality of a maximal independent set contained 
in Y. It is denoted rk(F). The closure of Y is defined as el(Y) := {x E A : 
rk(Y U x) = rk(Y)}. A set C C A is called a flat if C = cl(C). 

A set C C A is called a cocircuit if C n i? 7^ for all bases Bel and C is 
minimal with this property. Cocircuits of cardinality one are called coloops. An 
element that is not contained in any basis is called a loop. 

Now fix a basis Bel. An element b G B is called internally active in B 
if b = max(A \ cl(-B \b)), i. e. b is the maximal element of the unique cocircuit 
contained in (A\B) Lib. The set of internally active elements in B is denoted 1(B). 
An element x G A \ B is called externally active if x G c\{b <E B : b < x}, i. e. x is 
the maximal element of the unique circuit contained in BLix. The set of externally 
active elements with respect to B is denoted -E(-B)Q The Tutte polynomial 

T (AtB) (x,y):=^x^y lEm (1.2) 

BCI 

captures a lot of information about the matroid (A, B). One can show that it is 
independent of the order of the list A. 

In this paper, we mainly consider matroids that are realisable over some field K. 
Let X = (xi, . . . , xjv) be a list of vectors spanning some K- vector space W and let 
B(A) denote the set of bases for W that can be selected from X. One can easily 
see that (A, B(A)) is a matroid. The list A is called a realisation of this matroid 
and a matroid (^4, B) is called realisable if there is a list of vectors A and a bijection 
between A and X that induces a bijection between B and B(A). 

A standard reference for matroid theory is Oxley's book |37| . Survey papers on 
the Tutte polynomial are [H [23] . 

Throughout this paper, we use a running example, which we now introduce. 



Usually, combinatorialists use min instead of max in the definition of the activities. In the 
zonotopal algebra literature max is used. This has some notational advantages. 
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Example 1.1. 



Let X : = ( o i i ) = {xi,x 2 ,x 3 ). 



The set of bases that can be selected from X is M(X) = {(2:1,22), (xi,xa), (2:2, £3)}. 
The Tutte polynomial of the matroid (X,M(X)) is 

T(x^(X))(x,y) = x 2 + x + y. (1.4) 

1.3. Some commutative algebra. In this subsection, we define some commut- 
ative algebra terminology that is used in this paper. 

Definition 1.2 (A pairing between symmetric algebras). We define the following 
pairing: 

(...) :R[ Sl ,...,s r ] x E[ti, . . . , t r ] -> M (1.5) 

M*- ('(■&—■&)') ®- (L6) 

i. e. we let p act on / as a differential operator and take the constant part of the 
result. 

Remark 1.3. One can easily show that the definition of the pairing (■, ■) is inde- 
pendent of the choice of the bases for the symmetric algebras Sym([/) and Sym(V) 
as long as the bases are dual to each other. 

Definition 1.4. Let X C R[si,...,s r ] be a homogeneous ideal. Its kernel or 
Macaulay inverse system [551 [37J 13H is defined as 

kcrl := {/ G R[h,...,t r ] : (q, f) = for all q G 1}. (1.7) 
Remark 1.5. kerl can also be written as 

kerX:={fem[t 1 ,...,t r ]: P (J^,...,^-)f = 0} (1.8) 
where p runs over a set of generators for the ideal X. 

Remark 1.6. For a homogeneous ideal X C M[si, . . . , s r ] of finite codimension the 
Hilbert series of kerl and R[.si, . . . , s r ]/X are equal. For instance, this follows from 
H7J Theorem 5.4]. 

A graded vector space is a vector space V that decomposes into a direct sum 
V = (B;>o^4- A graded linear map f : V W preserves the grade, i. e. f(Vi) 
is contained in Wi. For a graded vector space, we define its Hilbert series as the 
formal power series Hilb(V,t) := ^ i>0 dim(Vi)i l . 

Note that a linear map / : V —> W induces an algebra homomorphism Sym(/) : 
Sym(V) -> Sym(W). 

1.4. Central zonotopal algebra. In this subsection, we define the Dahmen- 
Micchelli space V{X) and its dual V{X). The pair (D(X),V(X)), which is called 
the central pair of zonotopal spaces in [2S], is the origin of zonotopal algebra. 

A vector u G U naturally defines a polynomial p u G R[si,...,s r ] as follows: 
if u can be expressed in the basis (si, . . . , s r ) as u = AjSi, then we define 

Pu ■= Ya=i ^i s % e R [ s i> ■ ■ • , Sr]- For F C I, we define p y := l\ xeY Px- 
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Definition 1.7. Let X C U = W be a finite list of vectors that spans U . Then we 
define 

J(X) := ideal{p T : T C X cocircuit} C R[ Sl) . . . , Sr ] (1.9) 
andV(X) :=kcrJ(X) CR[t u ...,t r }. (1.10) 

2?(X) is called the central P-space or Dahmen-Micchelli space. 

It can be shown that T)(X) is the space spanned by the local pieces of the box 
spline and their partial derivatives. The box spline will be defined in the next 
section. The space V(X) was introduced in and in (TUJ it was shown that its 
dimension is |B(X)|. 

Definition 1.8. Let X C U = M. r be a finite list of vectors that spans U. Then, 
we define the central "P-space 

V{X) := span{p Y : Y C X,X \ Y has full rank} C R[s u ...,s r }. (1.11) 

Proposition 1.9 (|21j). Let X C U = M. r be a finite list of vectors that spans U . 
A basis for V(X) is given by 

B(X):^{Q B :BeB(X)} 1 (1.12) 

where Q B := P X \(bue(b))- 

The space V(X) can also be written as the kernel of an ideal. The following 
proposition appeared in [5] and earlier in a slightly different version in 

Proposition 1.10. Let X C U = M. r be a finite list of vectors that spans U. Then, 

V{X) = kerZ(X), (1.13) 

where 2(X) := ideal {p™^ : n G V \ {0}} C R[t u ...,t r ] (1.14) 

and m : V — > N assigns to i] G V the number of vectors in X that are not perpen- 
dicular to rj. 

Example 1.11. Let X be the list of vectors we defined in Example ll.il Then 

V{X) = kerideal{sis 2 , s\{si + s 2 ), s 2 (si + s 2 )} = spanjl, t\, t 2 }, (1-15) 
X(X) =ideal{i?,^,(ti - t 2 ) 2 } + R[ti,t 2 }> 3 = ideal{tl,tl,ht 2 }, (1.16) 
and V{X) = kerZ(X) = span{l, si, s 2 }. (1.17) 

Proposition 1.12 ([HI [3D]). Let X C U = M. r be a finite list of vectors that spans 
U . Then the spaces V{X) and T){X) are dual under the pairing (•, •), i. e. 

V(X)^V{Xy (1.18) 

f^(-J) (1-19) 

is an isomorphism. 

The preceding proposition implies that the Hilbert series of V{X) and T)(X) 
are equal. By Proposition 11.91 this Hilbert series is a matroid invariant and a 
specialisation of the Tutte polynomial. These facts are summarised in the following 
proposition. 
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Proposition 1.13. Let X C [/ = W be a list of vectors N vectors that spans U. 
Then 



mb(v(x),q) = m\b(v(x),q) = q N - r T (XMX)) (i^)= J2 v 



7 W-r-|_E(B)| 

BeR(X) 



11-20) 



Remark 1.14. Most of the results mentioned above also hold over other fields of 
characteristic zero or even over arbitrary fields. In this paper we work over the 
reals, because the construction in Section [3J uses distributions. Nevertheless, we 
believe that some of our results also hold over other fields. 

1.5. Organisation of the article. The remainder of this article is organised as 
follows. In Section [5] we provide some additional mathematical background, in 
particular on splines. 

In Section [3J we construct certain polynomials R B as convolutions of differences 
of multivariate splines. In Section [4] we show that 

B{X) :— {dct{B)R B : B 6 B(X)} (1.21) 

is a basis for T>(X) and we prove that this basis is dual to the basis B(X) for 
V(X). In Section[S]we discuss deletion-contraction and two short exact sequences. 
In Section [5] we introduce a new combinatorial structure called forward exchange 
matroid. This is an ordered matroid together with a subset B' of its set of bases with 
the so-called forward exchange property. In Section [7] we introduce the generalised 
P-space V{X,W) := span{Q B : B <G B'} and the generalised P-space V(X,W). We 
show that most of the results that we described in Subsection 11.41 and Section |4] still 
hold for these spaces if B' has the forward exchange property. For example, the two 
spaces are dual and a suitable subset of B(X) turns out to be a basis for T>(X, B'). 
Furthermore, T>(X,M') and V(X,M') have deletion-contraction decompositions that 
are related to the deletion-contraction reduction of the Tutte polynomial. 

In Section [8] we review the previously known zonotopal spaces and we show that 
they are special cases of our spaces T>(X,W) and V(X,W). 

Acknowledgements. The author would like to thank Olga Holtz and Martin 
Gotze for discussions and their helpful suggestions and Zhiqiang Xu for pointing 
the paper [43| out to him. 

2. Preliminaries 

In this section we provide some mathematical background. In Subsection 12.11 
we define some objects from discrete geometry. In Subsection 12.21 we review dis- 
tributions and in Subsection 12.31 we discuss box splines and multivariate splines. 
Subsection 12.41 contains information about previously known algorithms for the 
construction of bases for 2?-spaces. 

2.1. Cones and zonotopes. 

Definition 2.1. Let X = {x-y, . . . ,Xn) C U = M. r be a list of vectors. Then we 
define the zonotope Z{X) and the cone cone(X) by 

Z{X) := X i x i : < A ; < 1 1 and conc(X) := | ^ \ i x l : A; > o] . (2.1) 
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Zonotopes are closely connected to zonotopal algebra. For example, if X is 
totally unimodular, then vo\{Z(X)) = d\mV{X) = dim£>(AT) and the space V{X) 
can be obtained from the set of lattice points in the half-open zonotope via least 
map interpolation that is explained in Subsection 12.41 

2.2. Distributions. In this paper we are mainly interested in multivariate poly- 
nomials. However, in the construction in Section [3] more general objects appear 
as intermediate products. In this construction, we need "generalised polynomials" 
whose support is contained in a subspace. Furthermore, we use convolutions and 
the fact that convolutions and partial derivatives commute. Distributions have all 
of the desired properties. In this subsection we summarise important facts about 
distributions that we will need later on. For a detailed introduction to the subject, 
we refer the reader to Laurent Schwartz's book |39) . 

A distribution on a vector space U = M' (or an open subset of U) is a continuous 
linear functional that maps a test function to a real number. Test functions are 
compactly supported smooth functions [/->!. An important example of a dis- 
tribution is the delta distribution S x given by 5 x ((p) := tp(x). A locally integrable 
function /:{/—>• K defines a distribution Tf in the following way: 

T f (<p) := / f(u)<p(u)du. (2.2) 
Ju 

Recall that for two functions f,g:U — > R, the convolution is defined as 

f*9-= [ f(u)g(--u)du. (2.3) 
Ju 

This is well-defined only if / and g decay sufficiently rapidly at infinity in order 
for the integral to exist. The convolution of two distributions can also be defined 
under certain conditions. 

A distribution T vanishes on a set r C U if T(ip) = for all test functions whose 
support is contained in T. The support supp(T) of T is the complement of the 
maximal open set on which T vanishes. 

Let and T v be two distributions for which supp(S*^) n (K — supp(T r; )) is 
compact for any compact set K. Let cp : U ->lbea test function with support 
K. Then, we define the convolution 

(S 6 * T n )(<p) := S C (T„«M£ + r,))), (2.4) 

where a is a test function that is equal to one on a neighbourhood of supp(S'^) D 
(K — supp(T' I) )). When evaluating T v {a{C}^{£, + rj)), we think of ip(£_ + rj) as a 
function in rj and of £ as a fixed parameter. Then, T v (a(£,)(p(t; +rj)) is a function in 
£ with compact support that is contained in K — supp(T^). Note that the definition 
of («!?£ * Tr))(ip) is independent of the choice of the function a. The multiplication 
by a is necessary to ensure that T r/ (a(£)<p(£ + if)) as a function in £ has compact 
support. 

Note that the convolution of two distributions is a commutative operation and 
T * 6q = T. Let u G U . The partial derivative of a distribution T in direction 
u is defined by (D u T)(ip) := —T(D u ip). Convolutions of distributions have the 
same nice property with respect to partial derivatives as convolutions of functions. 
Namely, if T\ and are distributions on U and u € U, then 



D U {T X * T 2 ) = (A.T0 * T 2 = T x * (D U T 2 ). 



(2.5) 
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2.3. Splines. In this subsection we introduce multivariate splines and box splines 
as in [T?| Chapter 7]. Another good reference is |13) . 

Definition 2.2. Let X C U = W be a finite list of vectors. The multivariate 
spline (or truncated power) T x and the box spline Bx are distributions that are 
characterised by the formulae 

f(u)B x (u)du= I •••/ /( yXiXi ]dAi---dAjv (2.6) 



and 



f(u)T x {u)du = J ■J /(^AiXildAi-.-dAAr. (2.7) 



The multivariate spline is well-defined only if cone(JT) is pointed or equivalently 
if there is a functional <p <G V s. t. (p(x) > for all x £ X. It is of course always 
possible to multiply certain entries of the list X by — 1 s. t. this condition is satisfied. 
Note that in Definition 12.21 we do not require that X spans U in contrast to most 
of the rest of this paper. 

Bx and Tx can be identified with the functions 

B x( u ) = , \ f= vol7v-dim( sp a„(x)){2 6 [0; 1] N : X z = u} (2.8) 
\faet\XX 1 ) 

and T x {u) = ] = voljv_ dim ( spall (x)){2 G R> : = u}. (2.9) 

^/det(A X 1 ) 

It follows immediately from (|2.8p and (|2.9p that Bx is supported in the zonotope 
Z(X) and Tx is supported in the cone cone(X). For a basis C C J7, 

B c = and T c = (2.10) 

|det(C)| |det(C)| v ; 

Remark 2.3. The box spline can easily be obtained from the multivariate spline. 
Namely, 



Bx(x)= ^(-l)l s lT x (x-a s ), (2.11) 



sex 

where a s := Y^aes a - 

From now on, we will only consider Tx- We introduced the box spline only 
because of its importance in approximation theory. 

Theorem 2.4. Let X C U = W be a finite list of vectors that spans U . The cone 
cone(X) can be decomposed into finitely many cones Ci s. t. Tx restricted to each 
d is a homogeneous polynomial of degree N — r. 

Theorem 2.5. The space spanned by the local pieces of the multivariate spline Tx 
and their partial derivatives is equal to the Dahmen-Micchelli space T>(X) that was 
defined in Definition \1.7\ 

The multivariate spline can also be defined inductively by the convolution for- 
mula 



r ( x,x) - T x * T {x) = / T x {- - \x) dA 



(2.12) 
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FIGURE 1. The box spline and the multivariate spline defined by 
the list X in Example 11.11 The multivariate spline is calculated in 
Example 12.61 On the right, the dashed lines are level curves. 



using (|2.10p as a starting point. In particular, Tx = T Xl *- ■ -*T XN . Since D X T X = 5q, 
the convolution formula implies for Y C X that 

D Y T X = T X \ Y where D Y := JJ D x . (2.13) 

Example 2.6. We consider the same list X as in Example ll.il By (|2. 10|) . Tr XltX2 \ 
is the indicator function of M> . Then, by (|2.12| . we can deduce 

/>00 

T x (si,s 2 )= / Xr| («i - -M 2 - A)dA = mm(si,s 2 ). (2-14) 
<J o 

See Figure Q] for a graphic description of Tx ■ 
Example 2.7. Let X. t := (1, . . . , 1). Then, 

i times 

T Xl (s) = X M >0 (s) (2.15) 



and T Xi+1 (s) = / T Xi (s — A) dA = / - dA = — — — for a > 0. (2 



16) 



2.4. Previously known methods for constructing bases for P-spaces. Two 

other methods are known to construct a basis for T>{X). However, our algorithm 
has several advantages over the other two: it is canonical, i. e. it only depends on 
the order of the list X and it yields a basis that is dual to the known basis B(X) 
for the P-space. 

In Wolfgang Dahmen's construction [9] , polynomials are chosen as basis elements 
that are local pieces of certain multivariate splines. For certain choices of the 
parameters in his construction, it might yield the same basis as ours. 

The second construction uses the so-called least map interpolation that was 
introduced by Carl de Boor and Amos Ron [T3]. Given a finite set S C V, they 
construct a space of polynomials in 11(5) C Sym(V^) of dimension \S\ with certain 
nice properties. 
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Recall that U = R r , V denotes the dual space and a vector v G V defines a linear 
form p v G R[ii, . . . , t r ] = Sym(V). We define the exponential function as usual by 

e v ~Sr^ GK[[ti > ... J t r ]]SSym(lO*. (2.17) 

The least map J, maps a non-zero element of the ring of formal power series 
R[[ti, . . . , t r ]] to its homogeneous component of lowest degree that is non-zero. The 
least space of a finite set S C V is defined as 

U(S) := span{/j. : / G span{e 1 ' : v G S}} C R[i l7 . . . ,i r ]. (2.18) 

Let x £ U and G R. This defines a hyperplane 

H x , Cx := {veV: v(x) = cj. (2.19) 

If we fix a vector c G M. x , we obtain a hyperplane arrangement H(X, c) = {H x>Cx : 
x G X}. 

Every basis B C X determines a unique vertex 9b G V of the hyperplane 
arrangement H(X,c) that satisfies 6b{x) = c x for all x G -B. In matrix notation, 
Ob = B~ 1 cb, where cb denotes the restriction of c to M. B . If the vector c is 
sufficiently generic, then 9b ^ Ob' f° r distinct bases B and B' . In this case, 
the hyperplane arrangement H(X,c) is said to be in general position. For more 
information on hyperplane arrangements, see |40| . 

The following surprising theorem makes a connection between hyperplane ar- 
rangements and the space V(X). It generalises to other 2?-spaces (see [28l [29l [33] ) . 

Theorem 2.8 ([14)). Let X C U = W be a finite list of vectors that spans U. Let 
c G M. x be a vector s. t. the hyperplane arrangement H(X,c) is in general position 
and let S be the set of vertices ofH{X,c). Then 

V{X) = U{S). (2.20) 

De Boor and Ron give a method to select a basis from II(S') in [T5] (see also 
[T31 Chapter II] for a summary). Their construction depends on the choice of the 
vector c, an ordering of the bases and an ordering of N r while our construction only 
depends on the order on X. 

Example 2.9. This is a continuation of Example ll.il 

Let ci = c 2 = and c 3 = 1. The set of vertices of H(X, c) is 
S ={(0,0), (1,0), (0,1)}. Then, 

U(S) = span{/ ; : / G span{l, e l \ e' 2 }} 

= span{l, ti, £2}) since 1 = li, 

ti = (e* 1 - 1);, and t 2 = (e* 2 - 1);. 

3. Construction of basis elements 

In this section we construct a polynomial in R[.?i, . . . , s r ], given a finite list 
Z C U = R r and a basis BCZ. Later on we show that polynomials of this type 
form bases for various zonotopal 2?-spaces if one chooses suitable pairs (B, Z). The 
polynomial i?f is constructed as a convolution of differences of multivariate splines. 
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S 2 =R 2 




cone(P2) / 








Si 










bi = xi 






cone(— P2) 



FIGURE 2. The geometry of the construction of the polynomial 
R { x UX3} in Example H31 Note that N 2 = 0. 

Let Z C U be a finite list and let B = . . . , b r ) C Z be & basis. It is important 
that the basis is ordered and that this order is the order obtained by restricting the 
order on Z to B. For i S {0, . . . ,r}, we define Si = Sf := span{&i, . . . , bi}. Hence, 

{0} = Sjf C 5f C 5f C . . . C S B = [/ S R r (3.1) 

is a flag of subspaces. We define an orientation on each of the spaces Si by saying 
that (b%, . . . , bi) is a positive basis for Si. Now a basis D = (d±, . . . ,di) for S» is 
called positive if the map that sends b u to d v for 1 < v < i has positive determinant. 

Let u G Si \ Si— i. If (b\, . . . , bi—i, u) is a positive basis, we call u positive. 
Otherwise, we call u negative. We partition Z D (Si \ Si-i) as follows: 

P q B := {u e Z n (Si \ Si-!) : u positive} (3.2) 

and Nf := {u e Z (1 (Si \ Si-i) : u negative}. (3.3) 

We define 

Tf + := (-1)1^*1 • T Pi * T- Ni and Tf~ := (-l) lP ' 1 ■ T_ Pi * T Nt . (3.4) 
Note that Tf + is supported in cone(Pi, — JV») and that 

lf-(z) = (-l)^!^^). (3.5) 

Now define 

Rf := - Tf- and i?f = R B := R? * ■ ■ ■ * R B . (3.6) 

For an example of this construction see Example 14.41 and Figure^ In Corollarv l3.41 
we will see that the distribution R B can be identified with a homogeneous polyno- 
mial. 

Remark 3.1. A similar construction of certain quasi-polynomials in the discrete 
case is done in [XTl 13.6]. 

Remark 3.2. The construction of the polynomials R B may at first seem rather 
complicated in comparison with construction of the polynomials Qb that form 
bases of the P-spaces. 

Here are a few remarks to explain this construction: multivariate splines are very 
convenient because it is so easy to calculate their partial derivatives (cf. (|2.13l) ~). 
Taking differences of two splines in the definition of Rf ensures that R B is a 
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FIGURE 3. The setup in the proof of Proposition 13.31 

polynomial and not just piecewise polynomial. In fact, Rf * . . -*Rf is a "polynomial 
supported in S" for all i. 

We have to change the sign of some of the vectors before constructing the mul- 
tivariate spline T B+ to ensure that all the convolutions are well-defined. For ex- 
ample, the convolutions in (|3.6|) are well-dehned for the following reason: the sup- 
port of Rf * ■ ■ ■ * Rf is contained in Si . The support of Rf +1 is cone(Pi+i , — iVj+i ) U 
cone(— Pi+i, Ni+i). For every compact set K, the set 

{Si f)(K- (cone(P i+ i, -N i+1 ) U cone(P i+ i, -N i+1 )) (3.7) 

is compact. 

Proposition 3.3. The distribution R B is a local piece of the multivariate spline 
T?+*...*T r B +. 

Proof. Let c>0 and let 

t :=b 1 + -b 2 + ... + -^-oW-! + -^rb r . (3.8) 

c c c 

See Figure [3] for an example of this construction. The vector r is contained in 
cone(Z). By Theorem 12.41 there exists a subcone of cone(Z) that contains r s.t. 
Tz agrees with a polynomial p T ,z on this subcone. We claim that R B is equal to 
p T ,z- Note that 

Rz = {Ti + - Tf~) * • • • * (T, B + - T r B -) = (-l)W + 'Zi*j\ If i\+'Z ie j\P*\ Tzii 

JC[r] 

where 

Z J B = \J{P i ,-N i )u\J{-P i ,N i ). 

In order to prove our claim, it is sufficient to show that r is contained in conc(Z^) 
if and only if J = 0. The "if" part is clear. 

Let J be non-empty and let j* be the minimal element. For a <E K let <j> a : U — > R 
be the linear form that maps a vector x to 

r 

Y^(-l)x^) a i Xj , (3.9) 

3=3* 

where Xj denotes the coefficient of bj when x is written in the basis (pi,... ,b r ). 
We claim that for sufficiently large a, <j> a is non-negative on Z J B and <33 a (r) < 0. 
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By Farkas' Lemma (e.g. |38] Section 5.5]), this proves that r is not contained in 
cone(Z^). 

If x G Si n Z 3 B for i < j*, then obviously 4> a (x) = 0. If x G (Si \ n 
for i > j*, then A, ^ and \ v = for all v > i + 1. In addition, (-l) XJ «>Aj > 0, 
since all vectors in (Pi,— Ni) have a positive 6.; component when written in the 
basis (&!,.. . ,6 r ). Hence, = (-1) XJ (*V A, + o(a l ) = a* \X t \ + o(a i ). This is 

positive for sufficiently large a. For r, we obtain 

at?" a j * +1 a 3 * ( 1 \ 

— ^•"(^)- (3 ' 10) 

This is negative for sufficiently large c. Note that we fix a large a first and then we 
let c grow. □ 

Using Theorem [HH we can deduce the following Corollary. 

Corollary 3.4. The distribution R B can be identified with a homogeneous polyno- 
mial of degree \Z\ — r. 

Remark 3.5. The local pieces of the multivariate spline are uniquely determined by 
a certain equation (cf. |171 Theorems 9.5 and 9.7]). Taking into account Proposi- 
tion [33] this gives us a different method to calculate the polynomials R§ . 

The following theorem that is due to Zhiqiang Xu yields another formula for the 
polynomials R%. 

Theorem 3.6 (03] Theorem 3.1.]). Let X C U = W be a list of N vectors that 
spans U . Let c G WL X be a vector s. t. the hyperplane arrangement H(X, c) is in 
general position. For a basis B G M(X), let 9b G V denote the vertex ofH(X,c) 
corresponding to B (cf. Subsection \2A\) . Then 

Tx{u) = iN^y- B £) idet(^)iktr^- c ,) xconc(s)(u) - (3 - n) 

Note that the numerator p. lip is non-zero because c is chosen s. t. H.(X, c) is in 
general position. Using Proposition 13.31 one can deduce the following corollary. 



Corollary 3.7. Let Z CU = W be a list of n vectors that spans U. Let c and 9b 



as in Theorem 13.61 Then, the polynomial R B (u) is given by 



RBz{U) (»~r)l s , £ g B+) |det(50|n!^(^- C ,)' (3 ' 12) 

r £conc(i? / ) 

where r denotes the vector defined in (|3.8p and Z B+ denotes the reorientiation of 
the list Z s. t. all vectors are positive with respect to B, i. e. Z B+ = (Jl =1 (Pi, —Ni). 



4. A BASIS FOR THE DAHMEN-MlCCHELLI SPACE V(X) 

In this section we define a set B(X) and we show that this set is a basis for the 
central £>-space T>(X). Furthermore, we show that this basis is dual to the basis 
B(X) of the central "P-space V(X). Note that B is the equivalent of the letter B 
in the Cyrillic alphabet. 
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Definition 4.1 (Basis for T>(X)). Let X C U = W be a finite list of vectors that 
spans U . Recall that M(X) denotes the set of bases that can be selected from X 
and that E(B) denotes the set of externally active elements with respect to a basis 
B. We define 

B(X) := {dct(B)R^ E{B} : B € M(X)}. (4.1) 

Theorem 4.2. Let X C [/ = W be a finite list of vectors that spans U. Then 
B(X) is a basis for the central Dahmen-Micchelli space T>(X) and this basis is dual 
to the basis B(X) for the central V -space V(X). 

Remark 4.3. T>(X) and V(X) are independent of the order of the elements of X. 
The bases B(X) and B(X) both depend on that order. In Theorem 14.21 we assume 
that both bases are constructed using the same order. 

Example 4.4. This is a continuation of Example 11.11 See also Figure [5] The 
elements of B(X) are 



R Vu*,) = ^ _ T _ xi) # {T{x2 Xa) _ T{ _ X2 ,_ X3) ) = 8 2 , (4.3) 

and R [ ^ X3) = (T X2 - T_ X2 ) * (T (xi , Xa) - T ( _ X1> _ X3) ) = 8l . (4.4) 
The elements of B(X) are 

Q(x lt x 2 ) =P$ = 1, (4-5) 

Q(x l ,x 3 ) =Px 2 = h, (4.6) 

an d Q(x 2 ,x 3 ) = Pxi = h- (4.7) 
B(X) and B(X) are obviously dual bases. 



The proof of Theorem 14.21 is split into four lemmas. Recall that for a basis B = 
(b u . . . , b r ) we defined a flag of subspaces {0} = S o B CSfc...CS r B = £/^ W, 
where := span(6i, . . . , fej). 

Lemma 4.5 (Annihilation criterion). Let Z C U = W be a finite list of vectors 
and let B C Z be a basis. Let R% be the polynomial that is defined in (|3.6p . Let 
C C Z. Suppose there exists i £ [r] s. t. ZD (Sf \ SP_ 1 ) C C. 
Then, D C R§ = 0. 

Proof. Note that D a T_ a = —5q. Using (|2.13p . we obtain 

D c Rf = D c ((-l) lNil T Pi *T- Ni - (-l)\ p *\T- Pi *T Ni )) (4.8) 
= ^c\(s l \s l - 1 )(' 5 o-'5o) = 0. (4.9) 

This implies 

D C R B Z = Z? CV(SASi _ l} i?f * • ■ • * Rf_ x * * Rf +1 * • • • * R B = 0. (4.10) 

□ 

Lemma 4.6 (Inclusion). The polynomial Rx\e(b) * s contained in D(X \ E(B)) 
for all B e M(X). Since V(X \ E(B)) C V{X), this implies 

B(X) C (4.11) 
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Proof. Let det(B)R x * E , B , £ B(X) and let C C X \ E(B) be a cocircuit, i. e. C 
intersects all bases that can be selected from X \ E(B). We need to show that 
E>cRx\E(B) = 0. C can be written as C = X \ (H U E(B)) for some hyperplane 
H C U. 

Let i be minimal s.t. Si % H. Such an i must exist since S r = U. Even 
(Si \ n iJ = holds. This implies 

(X \ £(£)) n (Si \ Si-i) C X\(H U E(B)) = C. (4.12) 
By Lemma 1431 this implies DcRx\e(b) = ^. ^ 



The following lemma will be used only in the proof of Lemma 14.81 

Lemma 4.7. Let B,D £ B(X). Suppose that both bases are distinct but have the 
same number of externally active elements. 
Then there exists i £ [r] s. t. 

(X \ E(D)) n (SP \ SPJ a\(BU E(B)). (4.13) 

Proof. Let B = (bi, . . . , b r ) and D = (di, . . . , d r ). Suppose that the lemma is false. 
Then there exist vectors zi, . . . , z r s. t. 

Zi G (X \ E(D)) n (SP \SP_ 1 )n(BU E(B)). (4.14) 

Those vectors form a basis because z% £ SP \ SP_i- Since z% is not contained in 
E(D), Zi < di must hold. This implies E(D) C E(zi, . . . , z r ). On the other hand, 
E(zi, . . . , z r ) C E(B) since all Zj are contained in BU E(B). 

We have shown that E(D) C E(B). This is a contradiction since no finite set 
can be contained in a distinct set of the same cardinality. □ 

Lemma 4.8 (Duality). Let B,D £ M(X). Let Qb = Px\(bue(B)) € B(X) and let 
Rx\e(d) ^ e ^ e polynomial that is defined in (|3.6j) . Then 



5b, d denotes the Kronecker delta and we consider B and D to be equal if there 
exist 1 < i\ < . . . < i r < N s. t. B = (x^ , . . . , Xj r ) = D. 

Proof. By Corollary 13.41 Rx\e(d) 1S a homogeneous polynomial of degree N — 
r — \E(D)\. Thus, if |i?(-B)| ^ \E(D)\, then Q B and Rx\e(d) are homogeneous 
polynomials of different degrees and (Qb, Rx\eid)) = ^' 

Now suppose that B ^ D and both bases have the same number of externally 
active elements. In this case, the statement follows from Lemma [4.5l and Lemma [4.7l 

The only case that remains is B = D. Recall that R x \e(b) ~ R? * ■ ■ ■ * R? ■ 
Consider the ith factor Rf . The elements of (X \ E(B)) n (Si \ Sf_i) are used 
for the construction of Rf . Exactly one basis element is contained in this set: 
Recall that in Section[3]we defined a partition P, U N t = (X \ E(B)) n (Si \ Si-i). 
By construction, bi is positive, i. e. bi £ Pi. Now we apply the differential operator 
D(Pi\bi)\JNi to R% '■ 

D (PiUNi) \ bi ((-l)\ Ni \ -T Pi *T_ Ni - (-l)l^l . *TjvJ - (T 6i +T_ bj ). (4.16) 
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Now we can put things together. Note that X\(BU E(B)) = [f i=1 ((Pi \ h) U N t ). 
Hence, 

D x \(bue(b))Rb = (T bl +T_ 6l ) * • • • * (T br + T_ 6r ) = — . (4.17) 
This finishes the proof. □ 

Proof of Theorem 14.21 We know that V(X) and D(X) are dual via the pairing (•, •) 
and that B{X) is a basis for V{X). By Lemma B(X) and B(X) are dual to 
each other and by Lemma |4~6I B(X) is contained in D(X). Hence, B{X) is a basis 
for V{X). □ 

5. Deletion-contraction and exact sequences 

By Proposition 11.131 the Hilbert series of T)(X) and V{X) are equal and an 
evaluation of the Tutte polynomial. In particular, they satisfy a deletion-contraction 
identity that extends in a natural way to our algebraic setting. This is reflected by 
two dual short exact sequences. 

In this section we define deletion and contraction and we explain those two exact 
sequences. While the two sequences were known before, their duality has not yet 
been stated explicitly in the literature. 

Two important matroid operations are deletion and contraction. For realisations 
of matroids, they are defined as follows. Let X C U be a finite list of vectors and 
let x £ X. The deletion of x is the list X \ x. The contraction of x is the list X/x, 
which is defined to be the image of X \ x under the projection ir x :{/—>■ U jx. The 
space U and other terminology used here are defined in Subsection 11.11 

The space V(Xjx) is contained in the symmetric algebra Sym(£//x). If x = s r , 
then there is a natural isomorphism Sym([/ jx) = K[si, . . . , s r _i] that maps s,; to 
Si. This isomorphism depends on the choice of the basis (s±, . . . , s r ) for U. Under 
this identification, Sym(7r 2; ) is the map from R[si,...,s r ] to R[si, . . . , s r -i] that 
sends s r to zero and s\, . . . , s r -i to themselves. 

For T>(X/x), the situation is simpler: this space is contained in Sym((J7/x)*) = 
Sym(x°). This is a subspace of Sym(V). We denote the inclusion map by j x . If 
x = s r , then Sym((J7/x)*) is isomorphic to R[ti, . . . , t r -{\. This is a canonical 
isomorphism that is independent of the choice of the basis elements s±, . . . , s r _i. 

For a graded vector space 5*, we write S[l] for the vector space with the degree 
shifted up by one. 

Proposition 5.1 ([2])- Let X CU = K. r be a finite list of vectors that spans U and 
let x £ X be neither a loop nor a coloop. Then, the following sequence of graded 
vector spaces is exact: 

V(X \ x)[l] V{X) Sy "H x) V{X/x) -> 0. (5.1) 

Proposition 5.2 ( [16J ) . Let X C U = R r be a finite list of vectors that spans 
U and let x £ X be neither a loop nor a coloop. Then, the following sequence of 
graded vector spaces is exact: 

0^V(X/x) Hv(X)^V(X\x)[l}^0. (5.2) 

Note that (|5.2j) is a special case of (1.12) in [16] and it is exact by the results in 
that paper. 
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Remark 5.3. Proposition 15.11 and Proposition 15.21 are equivalent because of the 
duality of V(X) and V{X). 



Proof of Remark \5 .31 We only show that Proposition 15.21 implies Proposition 15.11 
The other implication is similar. 

Since dualisation of finite dimensional vector spaces is a contravariant exact 
functor, the following sequence is exact by Proposition 15.21 

o -> v(x \ x)* v{xy ( -H* v{x/ x y -> o. (5.3) 

By Proposition 11.121 V(X) is isomorphic to D{X)* via q n> (q, ■). Hence, it is 
sufficient to show that the following two diagrams commute: 

v(x\ x )^X?v(x\x)* v(x) 9 ^ {q 'i v(xy (5.4) 



►<?.•> 



and 



Sym(-7r x ) 



8 | -><?,->, 



> D(X)* —^V{X/x)* 

For the diagram on the left, we have to show that {p x q, •) = (<Z, fx ■ ) f° r a ll 
ge?(I\a;). This is easy. 

For the diagram on the right, we have to show that (SyTa(n x )q, ■) = (q,j x {-)) 
for all q £ V(X). If we choose a basis with s r = x this follows from the fact that 
^ r / = OforaU/GM[ti,...,t r _i]. □ 

6. Forward exchange matroids 

In this section we introduce forward exchange matroids. A forward exchange 
matroid is an ordered matroid together with a subset of its set of bases that satisfies 
a weak version of the basis exchange axiom (jl.ll) . 

The motivation for this definition is the following: the author noticed that most 
of the results in Subsection 1 1 . 41 and Section|4]hold in a far more general context. An 
important ingredient of the definitions of the spaces V(X) and T>{X) and their bases 
is the set of bases M(X) of the list X. These two spaces still have nice properties if 
we modify their definitions and use only a suitable subset B' of M(X ). It turned out 
that forward exchange matroids are the right axiomatisation of "suitable subset". 

Let (A,M) be an ordered matroid of rank r and let B = (pi, . . . ,b r ) £ B be 
an ordered basis. The flag (|3.1|l can be defined in combinatorial terms: for i £ 
{0, . . . , r}, we define Si = Sf :— c\{bi, . . . ,bi} C A. Hence, we obtain a flag of flats 

{x £ A: x loop} = 5,f C Sf C Sf C . . . C 5, B = A. (6.1) 

One can easily show that for a basis B £ B and i £ [r], the following statement 
holds: 

let x £ Sf \ S?_ v Then B' = (B \ b t ) U x is also in B. (6.2) 

Note that x £ Sf \ S?_ x satisfies x > bi if and only if x is externally active with 
respect to B. This motivates the name of the following definition. 

Definition 6.1 (Forward exchange property). Let (A, B) be a matroid and let 
B' C B. We say that the set of bases B' has the forward exchange property if the 
following holds for all bases B £ B' and all i £ [r]: 

let x £ Sf\ (SfL ! U E{B)). Then B' = (B \ b t ) U x is also in B'. (6.3) 



18 



MATTHIAS LENZ 



Remark 6.2. Note that Sf = Sf holds for all j > i. If x is the ith vector in B', 
this equality holds for all j, i. e. B and B' define the same flag. 

Definition 6.3 (Forward exchange matroid). A triple (A,M,M') is called a forward 
exchange matroid if ( A, B) is a matroid and B' is a subset of the set of bases B with 
the forward exchange property. 

Remark 6.4. In this paper, we mainly consider realisations of forward exchange 
matroids, i. e. pairs (X,W) where X is a list of vectors and B' C M(X) is a set of 
bases with the forward exchange property. 

Definition 6.5 (Tutte polynomial for forward exchange matroids). Let (A, B,B') 
be a forward exchange matroid. 

Then, we define its Tutte polynomial to be 

T {AtMiW) (x,y):= T£ x \W\ y \*W\ t ( 6 .4) 

BGB' 

where 1(B) and E(B) denote the sets of internally and externally active elements 
with respect to B in the matroid (A,M). 

Remark 6.6. It would be interesting to clarify the relationship between forward- 
exchange matroids and other set systems studied in combinatorics such as greedoids 

IH El]. 

7. Generalised £>-spaces and "P-spaces 

Earlier, we considered the spaces T>(X) and V(X) for a given list of vectors X. 
The construction of those spaces relied mainly on the matroidal properties of the 
list X, namely on the sets of bases and cocircuits. 

Motivated by questions in approximation theory, various authors generalised 
those constructions. Given a list X and a subset B' of its set of bases M(X), one can 
define a set T*(X,W) as the kernel of the ideal generated by the B'-cocircuits (i. e. 
sets that intersect all bases in B'). Under certain conditions, dimP(X, B') = |B'| 
still holds. In this section we show that if the set B' has the forward exchange 
property, this equality holds and there is a canonical dual space V(X,M'). Both, 
the generalised 2?-spaces and the generalised 7^-spaces satisfy deletion-contraction 
identities as in Section [5] and there are canonical bases for both spaces that are 
dual. 

7.1. Definitions and Main Result. 

Definition 7.1 (generalised 2?-spaces). Let X C U = R r be a finite list of vectors 
that spans U and let B' be an arbitrary subset of its set of bases B(A). A set CCI 
is called a B' -cocircuit if C intersects every basis in B' and C is inclusion-minimal 
with this property. 

The generalised T>-space defined by X and B' is 

V{X, B') := {/ : D c f = for all B'-cocircuits C} = kcr J(X, B'), (7.1) 

where J(X, B') := ideal{p c : C C X is a B'-cocircuit}. (7.2) 

Proposition 7.2 ( |14( Theorem 6.6]). Let X C U = R r be a finite list of vectors 
that spans U and let B' be an arbitrary subset of its set of bases M(X). Then 

dimV(X,W) > |B'|. (7.3) 
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Definition 7.3 (generalised "P-spaces). Let X C U = M. r be a finite list of vectors 
that spans U and let B' be an arbitrary subset of its set of bases M(X). Then we 
define 

B(X, B') := {Q B ■ B G B'} = fe\ (BU £(B)) : B G M'} (7.4) 

and P(X,B') := spanB(X,B'). (7.5) 

We call V(X,M') the generalised V -space defined by X and B'. 

Remark 7.4. The set B(X,M') is a basis for V(X, B'). By definition, it is spanning 
and it is linearly independent because it is a subset of B(X). 

If the set B' has the forward exchange property, the spaces T)(X, B') and V(X, B') 
have many nice properties. Here is the Main Theorem of this section. 

Theorem 7.5. Let X C U = K. r fee a finite list of vectors that spans U and let 
B' C B(X) fee a se£ o/ feases with the forward exchange property. 

Then, the generalised V-space D(X,M') and the generalised V -space V(X,M') are 
dual via the pairing (•,•). In addition, 

B(X,W) := {det(B)R^ E{B) : B G B'} (7.6) 

is a basis for D(X,M') and this basis is dual to the basis B(X,M') for V{X,W). 

Corollary 7.6. Let X C U = W be a list of N vectors that spans U and let 
B' C M(X) be a set of bases with the forward exchange property. Then 

mib(v(x,W),q) = mib(v(x,m'),q) = ^ ^-r-^m = q ^ Tmx)<w) {i, -). 

BEM' " 

Here are two examples that help to understand generalised D-spaces, generalised 
■P-spaces, and Theorem 17.51 

Example 7.7. Let X — [ex, e2, a, b) C M 3 where ei,e2,e3 denote the unit 
vectors and a = (a, ft, 7) and fe are generic. In particular, a, f3, 7 ^ 0. Let 

B' := {(eie 2 e 3 ), (eie 2 a), (eie 2 fe), (eie 3 o), (eie 3 6), (e 2 e 3 a), (e 2 e 3 fe)} C M(X). 

The reader is invited to check that B' has the forward exchange property. The 
Tutte polynomial is T(x,M(x),m')(x,y) = 3a; + 3y + y 2 and the B'-cocircuits are 
{eie 2 , eie 3 , e 2 e 3 , eiafe, e 2 afe, e 3 afe}. Hence, 

T>(X,W) = kcrideal{sis 2 , S1S3, s 2 s 3 , sip ab , s 2 p a b, s 3 p ab } 

= span{l, ti, t 2 , t 3 , t\, t\, tl}, 

B(X,W) = {1, p e3 ,Pe 3 a, Pe 2 , Pe 2 a, P ei , Pe ia }, 

V(X,W) = span{l, si, s 2 , s 3 , Si(asi + /3s 2 +7^3), s 2 (asi + f3s 2 +7 S 3), 

s 3 (as 1 + /3s 2 +^s 3 )}, and 

B(XX) = ll,t 3 , ^,t 2 , 4,*i, — V 

Example 7.8. Let iV > 3 be an integer. Let = (xi, . . . , xjv) be a list of 
vectors in general position in R 2 with x\ = e\, x 2 =■ e 2 , and x 3 = e\ + e 2 . In 
addition, we suppose that the second coordinate of all vectors Xi (i > 3) is one. Let 
B' := {(xi,Xi) : i G {2, . . . , iV}}U{(a; 2 , x 3 )}. Note that B' is totally unimodular, i. e. 
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all elements have determinant 1 or —1 and W has the forward exchange property. 
Then, 

T>(X N ,W) = kerideal{p XlX2 ,p XlX3 ,p X2 ... XN } = span{l,ti,* 2 , ■ ■ - ,^~ 2 }, 

f t t N ~ 2 1 

^(Xjv, BO = | 1, fr, t 2 , |, . . . , (jV 2 _ I , and (7.7) 

B(X N ,W) = {l, Pxi }U{p X2 ... Xi :ie{2,...,N-l}}. (7.8) 



Now we embark on the proof of Theorem 1 7. 5 1 We start with the following simple 
lemma. 

Lemma 7.9 (Inclusion). Let X C U = R r &e a finite list of vectors that spans U 
and let B' C B(X) fee a se£ o/ 6ases with the forward exchange property. Then 

B(X,W) <ZV(X,W). (7.9) 

Proof. Let Bel' and let det(B)i?J^ B ^ G B(X,M') be the corresponding basis 
element. Let C C X be a B'-cocircuit, i. e. an inclusion-minimal subset of X that 
intersects every basis in B'. 

Let B = (bi, b r ). If there exists an i s. t. (S*f \ S* 3 ^) f) (X \ E(B)) C C, 
we are done by Lemma 14.51 Now suppose that this is not the case, i. e. for every 
i S [r], there is a e (5f \Sf_ 1 )n(X\ (E(B) U C)). Then we define a sequence 
of bases B , ■ ■ ■ , B r by 

B :=B and := (^_i \ 6,-) U Zi for i G [r]. (7.10) 

The lists _Bj are indeed bases and even though in general, they might define different 
flags, they satisfy 

(Sf - 1 \ n (X \ B^-O) = (Sf \ Sf_ x ) n (X \ E(B)), (7.11) 

because span(6i, . . . ,bi) = span(zi, . . . , Zi) for all i g [r]. Hence, Bi G B' implies 
G B' because B' has the forward exchange property. In particular, B r = 
[z\, ... , z r ) G B'. By construction, B r DC = 0. This is a contradiction. □ 

Definition 7.10. Let (A,M) be a matroid and let B' C B. Let x e A. B' can be 
partitioned as B' = lu U Bu, where 

M[ x := {B G 1' : x <£ B} denotes the deletion of x and (7.12) 

Bj^ := {B G 1' : x G 5} the restriction to x. (7.13) 

If we are given a list of vectors X C U and a set of bases B' C B(X), we can 
also define the contraction W, x . Recall that ir x ; U U / x denotes the canonical 
projection. Then, we define 

W /x :={n x (B\x):xEBeM\ x }. (7.14) 

Remark 7.11. For technical reasons, it is helpful to distinguish the contraction M'^ x 
and the restriction Bj^, although there is a canonical bijection between both sets. 

We now introduce the concept of placibility. This is a condition on a set of bases 
B' which implies equality in (|7.3p . 

Definition 7.12 ([IB], see also [33]). Let (A, B) be a matroid and let B' C B be a 

non-empty set of bases. 
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(i) We call an element x & A placeable in B' if for each Bel', there exists an 
element be B such that (B \ b) U x G B'. 

(ii) We say that B' is placible if one of the following two conditions holds: 

(a) B' is a singleton or 

(b) there exists x G A s. t. x is placeable in B' and both, B'^ and Mi x are 
non-empty and placible. 

Proposition 7.13 ([16]). Let X C U = M. r be a finite list of vectors that spans 
U and let B' be an arbitrary subset of its set of bases M(X). J/B' is placible, then 
dimP(X,B') = |B'|. 

Lemma 7.14. Let (A,M,W) be a forward exchange matroid. ThenM' is placible. 

Proof. If |B'| = 1, then B' is placible by definition. Now let |B'| > 2. Let x be the 
minimal element in A s. t. both, B(, and B'^ are non-empty. Such an element must 
exist if |B'| > 2. 

We show now that x is placeable in B'. Let B = (bi, . . . , b r ) be a basis in B' 
and let i G [r] s. t. x G Sf \ Sf_ 1 . We claim that x < bi. Suppose not. Because 
of the minimality of x, this implies that bi, . . . , bi are contained in all bases in B'. 
Since x G span(6i, . . . , bi), this implies that x is not contained in any basis. This is 
a contradiction because we assumed that B'^, is non-empty. 

Now we have established that x < hi. This implies that x is not externally 
active. Hence, because of the forward exchange property, (B \ bi) U x G B', i. e. x 
is placeable in BeB'. 

It remains to be shown that Bt, and B'^ are both placible. By induction, it is 
sufficient to show that both sets have the forward exchange property. For B( ^ , this 
is clear. For Bj , this follows from the following fact: by the choice of x, all a G A 
that satisfy a < x are either contained in all bases or in no basis in Bf. □ 

Proof of Theorem\TM By Lemma EH B(X,M') C T>(X,W). By Proposition EH 
and by Lemma EH dim£>(X,B') = |B'| = \B(X,W)\. Linear independence of 
B(X,M') is clear because it is a subset of B(X). For the same reason, the duality 
with B(X, B') C B(X) follows from Lemma EU □ 

Remark 7.15. The correspondence between T>(X) and the set of vertices S of a 
hyperplane arrangement l-L{X,c) in general position that is stated in Theorem 12.81 
generalises in a straightforward way to a correspondence between V(X,M') and the 
subset of S that is defined by B'. 

7.2. Deletion-contraction and exact sequences. In this subsection, we show 
that the results in Section [S] about deletion-contraction and exact sequences natur- 
ally extend to generalised 2?-spaces and "P-spaces. We use the same terminology as 
in that section. 

Recall that for a graded vector space S, we write S[l] for the vector space with 
the degree shifted up by one. 

Proposition 7.16. Let X C U = K r be a finite list of vectors that spans U and 
let B' C M(X) be a set of bases with the forward exchange property. Let x be the 
minimal element of X that is neither a loop nor a coloop. Then, the following 
sequence of graded vector spaces is exact: 

-> T>(X/x,m' /x ) h T>(X,W) ^ V{X \ x,M[ x )[1] -> 0. (7.15) 
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Proof. Let B e B' be a basis that does not contain x. Because of the minimality, x 
is not externally active with respect to B. This implies D x Rx\e(b) = ^X\(E(B)Uxy 
Hence, D x : {dct(B)R B g E(X,W) : x £ B} -> S(X \ a;,B' ) is a bijection and 
consequently, Z) x maps T>(X,W) surjectively to V(X \ x,Mi). For the rest of the 
proof, we refer the reader to [16], in particular to the explanations following (1.12) 
and to Theorem 2.16. □ 

Proposition 7.17. Let X C U = M. r be a finite list of vectors that spans U and 
let B' C M(X) be a set of bases with the forward exchange property. Let x be the 
minimal element of X that is neither a loop nor a coloop. Then, the following 
sequence of graded vector spaces is exact: 

o -> v(x \ x,m[ x )[i] ^ p(x,W) Sy ^ x) v(x/x,m' /x ) ->. o. (7.16) 

Proof. One can easily check this for the bases of the "P-spaces. Alternatively, it 
can be deduced from Proposition 17.161 using a duality argument as in the proof of 
Remark 15.31 □ 



Remark 7.18. The exact sequences in this section require x to be minimal in contrast 
to the ones Section[Sl where x can be any element that is neither a loop nor a coloop. 
This reflects the fact that matroids have an (unordered) ground set, while forward 
exchange matroids have an (ordered) ground list. 

Remark 7.19. One could replace V{X/x,W /x ) by T>(X,W.) in Proposition OH 
The analogous replacement in Proposition 17 17l would be problematic. The reasons 
for that are explained in Section [3] 

7.3. VlX, B') as the kernel of a power ideal. By now, we have seen that most of 
the results about T>{X) and V(X) that we stated earlier also hold for the generalised 
2?-spaces and P-spaces. The only thing that is missing is a power ideal X(X,M') 
s. t. V(X,M') = kerI(X, B'). Unfortunately, such an ideal does not always exist. 

In this section, we describe the natural candidate for this power ideal and we 
give an example where its kernel is equal to V{X, B') and one where it is not. 

Definition 7.20. Let X C U = W be a finite list of vectors that spans U and let 
B' be an arbitrary subset of its set of bases M(X). Recall that V := U* . We define 
a function K : V — > N by 

k(t/) := m^\X\(BU E(B)Uif)\ (7.17) 
and 1(X, B') := ideal{p«W +1 : r, G V \ {0}}. (7.18) 

Lemma 7.21. Let X C U = W be a finite list of vectors that spans U and let B' 
be an arbitrary subset of its set of bases M(X). Then 

V(X, 1') C kcrI(X, B'). (7.19) 

Proof. It is sufficient to show that all elements of the basis B(X, B') are contained 
in kcrX(X,B'). Let Bel' and let n e V \ {0}. Then 

D^Px\(byje(b)) = Pixn-n^MBuEiB^D^^PxXiBuEiB^v") = (7.20) 

The first equality follows from Leibniz's law. The second follows from the fact that 
by definition, k(t]) >\X\(BU E{B) U □ 
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Remark 7.22. If one examines the proof of Lemma T7. 2 11 one immediately sees that 
I(X, B) is the only power ideal for which V(X, B') = keiI(X, B') can possibly hold. 

Remark 7.23. In some cases, V(X, B') and kerX(X, B') are equal (see Example l7.25ll . 
In other cases however, V(X, B') is not even closed under differentiation (see Ex- 
ample o. 

Remark 17.231 naturally leads to the following question. 

Question 7.24. Is there a simple criterion to decide whether V{X 1 B') is closed 
under differentiation or if V(X,3') = kerJ(X,B') holds? 

Example 7.25. This is a continuation of Example 17.71 Recall that we considered 
the list X = (ei, e2, e3, a, b) C R 3 where a and b are generic vectors and a = (a, /3, 7) 
with a, (3, 7 ^ 0. The set of bases is 

B' = {(eie 2 e 3 ), (eie 2 a), (eie 2 6), (eie 3 a), (eie 3 6), (e 2 e 3 a), (e 2 e 3 6)} C1(I). 

In order to calculate the function k, we first determine the inclusion-maximal lists 
in {X \ (B U E{B)) : B e B}. Those are (eia), (e 2 a), and (e 3 a). We can deduce 
that k(t]) is one if rj £ a° and two otherwise. We obtain 

1(X,W) = idcal{p( Qi _ ( g !0) ,p( 0!/3 _ 7) ,p( Qi0i _ 7) } + M[si,s 2 ,s 3 ]> 3 and 

V(X,W) =kerI(X,B') = span{l, s 1; s 2 , s 3 , Sl (a Sl + £s 2 + 7 s 3 ), 

s 2 (asi + /3s 2 + 7S3), s 3 (asi + /3s 2 + 7s 3 )}. 

The degree two component of M[s 1; ,s 2 , s 3 ]/I(X, B') is three-dimensional. This im- 
plies that V(X,W) = keiT(X,W). 

8. Comparison with previously known zonotopal spaces 

In this section we review the definitions of various zonotopal spaces that have 
been studied previously by other authors. It turns out that they are all special 
cases of the generalised P-spaces and P-spaces that we introduced in Section [7] 
The most prominent examples are of course the central spaces T>(X) and V(X) 
that we obtain if we choose B' = M(X). 

Let A = (oi,...,o n ) C l r be a list of vectors that spans W and let X = 
(xi, . . . , xn) C W\ where N > n and a.; = Xi for i £ [n]. In |28l HH1 [S3, the spaces 
T>(X,W) and V(X,B') are studied for certain sets of bases B' C M(X). Here are 
the definitions of those sets of bases. 

Definition 8.1 (Internal and external bases [21] )• Let A C W be a list of vectors 
that spans R r and let Bo = (61, . . . , b r ) C W be an arbitrary basis for K r that is 
not necessarily contained in B(^4). Let X = (A, Bq) and let 

ex: {I C A: I linearly independent} -> M(X) (8.1) 

be the function that maps an independent set in A to its greedy extension. This 
means that given an independent set / C A, the vectors b\,.. . ,b r are added suc- 
cessively to / unless the resulting set would be linearly dependent. 

Then we define the set of external bases B + (A, Bq) and the set of internal bases 
M-{A) by 

M + (A, B ) := {B e M(X) : B = ex(I) for some I <ZA independent} (8.2) 
and B_(^4) := {B G M(A) : B contains no internally active elements}. (8.3) 
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The lattice of flats C{A) of the matroid (A, B(A)) is the set {C C A : cl(C) = C} 
ordered by inclusion. An upper set J C is an upward closed set, i. e. C\ C 

C2 G J implies Ci G J. 

Definition 8.2 (Semi-internal and semi-external bases [29]). We use the same 
terminology as in Definition 18. f I In addition, we fix an upper set J in the lattice 
of flats C(A) of the matroid (A,~B(A)). For the semi-internal space, we fix an 
independent set Iq C A whose elements are maximal in A. 

Then we define the set of semi- external bases M + (A, B , J) and the set of semi- 
internal bases M_(A,I ) by 

M + (A, B , J) := {B G M(X) : B = cx(I) for some I C A independent 

and cl(J) e J} and 

B_(A, Iq) := {B G M(A) : B (~l Iq contains no internally active elements}. 

Definition 8.3 (Generalised external bases [33]). Let A = (01, . . . , a„) C K r be a 
list of vectors. Let k : C(A) —> {0, 1,2,.. .} be an increasing function, i. e. C\ C C2 
implies k(C\) < n(C2)- 

Let X = (j4, Y), where Y = (j/Xj 2/2, • ■ • , y K (A)+r) is a list of generic vectors, i. e. 
if yi is in the span of Z C X \ j/j, then span(Z) = span(X). 

Then we define 

B K (A,F) :={BGB(X) :Snr C (y 1 ,...,y K(cl(AnB))+ | Bny |)}. (8.5) 

Remark 8.4. The spaces V{X,W) and X»(X,B') are equal to 

• the external spaces V+(X) and V + {X) in if B' is the set of external 
bases; 

• the semi-external spaces V + (X, J) and V + (X, J) in [5S] if B' is the set of 
semi-external bases; 

• the generalised external spaces V K {X) and T> K {X) in [33J if B' is the set of 
generalised external bases. For V K (X), we need to assume in addition that 
k is incremental, i. e. for two flats C\ C C2, k.(C<x) — k{C\) < dim(C2) — 
dim(Ci). 

Furthermore, the space T>(X,M') is equal to the (semi-) internal space 2?_(X) resp. 
D-(X,I ) in [5H1IM] if B' is the set of (semi-)internal bases. 

Remark 8.5. The (semi-)internal spaces V{X,M-(X)) and V(X, B_(X, J )) are in 
general different from the spaces V-{X) and V-(X,Io) in [28] [55], but they have 
the same Hilbert series. 

Remark 8.6. The theorems about duality of certain 'P-spaces and 2?-spaces in |28l 
[551 [35] are all special cases of Theorem 17.51 This is a consequence of Lemma 18.71 
below. 

Lemma 8.7. The sets of bases defined in Definitions \8A[ 18.21 and \H.'S\ all have the 
forward exchange property. 

Proof. We use the following notation throughout the proof: B = (bi, . . . , b r ) is a 
basis and x G (Sf \S ? i3 _ 1 )n(X\i;(.B)) for some i G [r]. In addition, B' := (BUx)\bi. 
Since x is not externally active, x < bi holds. We may even assume x < bi because 
if equality occurs, nothing needs to be shown. 

Internal and external bases are special cases of semi-internal and semi-external 
bases so we do not consider them separately. 
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We start with the semi-external bases. Let B G M+(A, Bo, J), i. e. B is the 
greedy extension of an independent set I C A. Recall that x < b%. Hence, x G A 
because if x was in Bq, the greedy extension of I would contain x instead of 6j. 
Now one can easily check that B' n A is independent and that ex(B' D A) = B' . 
This is equivalent to B' G M + (A, B ,J). 

Now we consider the semi-internal bases. Let B G M-(A, Jo). By construction, 
the fundamental cocircuits of B and bi resp. x are equal. As x < hi, inactivity of 
hi implies inactivity of x. Hence, B 1 G M-(A,Iq). 

Last, we consider the generalised external bases. Let B G M K (A,Y). If bi G A, 
then K (c\(AnB)) + \B nY\ = n(cl(Ar)B')) + \B' DY\. This implies B' G M K (A, Y). 
If 6j G Y , then j = r must hold because the vectors in Y" are generic and we are 
supposing bj ^ x. Replacing b r by x reduces the index of the maximal element in 
Y that is permitted in B by at most one since k is non-decreasing on C{A). Since 
we remove the maximal element of the basis B, this causes no problems. □ 

Remark 8.8. Various "P-spaces without a dual £>-space have been studied by several 
other authors e.g. [3J EH ES3 EH- The approach in d [32] is slightly different 
from ours. The authors of these two papers start with a list of vectors X and an 
integer k. In our construction, it is eventually necessary to add additional elements 
to the list X in order to obtain arbitrarily large "P-spaces. In their construction, it 
is sufficient to let the integer parameter k grow while keeping the list X fixed. 

Their construction of a basis for the P-space takes into account the internal 
activity of the bases in B(X). Every element of M(X) with internally active elements 
may define multiple elements of the basis for the P-space. 

Example 8.9. Example 17.81 fits into the framework of [33] resp. Definition 18.31 if 
we choose A = (xi), Y = (x2, ■ ■ ■ ,xn) and n as follows: «(span(ei)) := N — 2 and 
on all other one-dimensional flats C, n{C) := 0. In this case, our space T>{X,W) 
is the same as the space T> K in [33]. The P-spaces are different because k is not 
incremental. 
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